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ATMOSPHERIC  PHOTO  IONIZATION  OF  ARGON  WITH  A  SHORT  PULSE  LASER 


I.  Introduction: 

One  of  the  goals  of  the  NRL  laser/atmosphere  program  is  to  use  a  short  pulse  laser  to  produce  an 
atmospheric  breakdown  at  a  controlled  distance  from  the  laser.  The  currently  favored  way  of 
doing  this  is  to  use  the  dispersive  properties  of  the  atmosphere  to  compress  a  pulse  longitudinally 
until  it  gets  to  the  irradiance  where  self  focusing  can  compress  it  in  the  transverse  plane  [1],  The 
compressed  pulse  reaches  such  an  irradiance  that  it  generates  breakdown  in  the  atmosphere. 
Although  final  laser  parameters  are  not  specified  at  this  point,  as  a  rule  of  thumb,  we  think  of  a  1 
ps  pulse,  that  compresses  to  100  fs  as  it  propagates,  and  we  think  of  the  irradiance  after  self  focus 
as  being  about  lO'^  -  lO'^  W/cm^.  Clearly  to  control  the  distance  of  breakdown,  one  must 
control  the  bandwidth  of  the  laser;  as  the  spread  in  group  velocity  increases,  the  focal  point 
moves  closer  to  the  laser. 

However  the  chirped,  short  pulse  lasers  typically  radiate  only  about  one  Joule  per  pulse.  Taking 
an  ionization  energy  of  15eV,  this  can  generate  at  most  4xl0'^  electron  ion  pairs  if  no  laser 
energy  is  lost  to  other  processes.  Furthermore,  in  the  atmosphere,  these  electrons  last  only  a  few 
ns  due  to  three  body  attachment  to  oxygen.  For  certain  applications,  this  may  be  sufficient. 

Other  applications  might  require  more  electrons  and/or  longer  life.  For  these  latter  applications, 
the  photo  ionized  plasma  could  serve  as  a  seed  plasma.  This  seed  plasma  could  be  ionized  by  for 
instance  a  focused,  high  power  CO2  laser.  As  long  as  the  seed  electrons  are  present  in  the 
atmosphere,  a  CO2  laser  is  ideal  (right  frequency,  high  power,  high  efficiency)  for  avalanche 
breakdown  of  the  atmosphere.  However  the  frequency  of  the  CO2  laser  is  far  too  low  to  generate 
the  required  seed  ionization  by  multi-photon  seed  ionization  (it  would  be  more  than  a  100  photon 
process),  thereby  opening  an  opportunity  for  seed  ionization  by  a  much  higher  frequency  laser. 

The  pulse  time  of  this  high  frequency  laser  is  too  short  for  avalanche  ionization,  and  the 
irradiance  is  too  low  for  tunneling  ionization.  Hence  one  is  left  with  multi  photon  ionization. 

The  question  is  how  to  pick  the  laser  parameters  so  as  to  optimize  this  process.  Before  we  begin 
we  briefly  discuss  photo  ionization.  If  a  photon  has  energy  above  the  ionization  energy,  photo 
ionization  usually  proceeds  very  quickly  for  a  high  power  laser.  However,  if  its  energy  is  below 
the  ionization  energy,  two  or  more  photons  may  combine  coherently,  or  incoherently  to  ionize  the 
atom.  Multi-photon  ionization  has  been  a  well  known  process  for  a  long  time  [2-5].  Numerous 
experiments  have  studied  it  [6].  If  the  ionization  energy  is  Ej  and  the  photon  energy  is  Eph,  one 
would  naturally  think  that  minimum  number  of  photons  required  would  by  Ej/Eph.  However 
because  of  the  presence  of  excited  states  below  the  ionization  threshold,  it  is  also  possible  that 
that  fewer  photons  could  be  required  because  the  ionization  might  proceed  via  a  cascade  process 
through  one  or  more  of  these  excited  states.  Such  a  process  has  been  called  resonantly  enhanced 
multi  photon  ionization  (REMPI).  Recently  this  author  has  argued,  using  experimental  data  and 
simple  theoretical  concepts  [7],  that  for  short  pulse,  high  power  lasers,  REMPI  could  occur  quite 
frequently,  indicating  that  the  number  of  photons  required  could  be  more  accurately  estimated  as 
something  like  (0.75-0.9)  Ej/Eph.  This  memo  suggests  one  particular  strategy  to  optimize 
REMPI  in  the  atmosphere. 

It  suggests  that  ionizing  the  argon,  rather  than  the  oxygen  or  nitrogen  could  be  the  easiest  way  to 
generate  the  ionization.  The  atmosphere  is  about  78%  N2,  about  21%  O2,  and  about  1%  Ar. 
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Thus  there  is  plenty  of  argon  to  produce  a  dense  plasma.  The  ionization  energy  of  O2  is  12.1  eV, 
that  of  N2  is  15.6  eV,  and  that  of  Ar  is  15.8  eV,  so  one  could  naturally  question  why  the  argon 
would  ever  be  the  easiest  component  to  ionize.  The  possible  answer  may  be  that  argon  is  an 
atom,  whereas  the  other  two  are  molecules.  Thus  the  intermediate  electronic  excitations  of  argon 
are  isolated  states,  whereas  the  excited  electronic  states  of  the  molecule  occur  in  broad  bands 
because  of  the  presence  of  so  many  nearby  vibrational  and  rotational  states.  As  we  will  discuss 
shortly  (and  as  is  well  known),  the  excitation  of  isolates  states  is  usually  a  much  stronger  process 
than  excitation  states  in  a  band.  Thus,  if  the  laser  frequency  is  picked  properly,  one  might  be 
able  to  ionize  the  argon  component  of  the  air  by  exciting  through  this  intermediate  state, 
furthermore,  this  intermediate  state  may  be  easier  to  access  than  the  analogous  intermediate 
molecular  state.  Others  have  also  suggested  this  same  approach,  for  other  applications,  using 
both  atmospheric  argon  [8]  and  xenon  [9]. 

Let’s  now  scope  out  the  problem.  We  want  the  photon  energy  to  be  as  high  as  possible,  so  as  few 
photons  as  possible  are  required  for  the  multi  photon  processes.  However  the  photon  must  also 
propagate  freely  in  the  atmosphere.  The  shortest  wavelength  photon  that  propagates  essentially 
loss  free  in  the  atmosphere  has  a  wavelength  of  about  3750  Angstroms,  or  about  3.3  eV.  The 
lowest  energy  electronic  excitation  of  argon  (designated  4s[l  V2]°)  has  an  energy  of  about  1 1.65 
eV  above  the  ground  state  (designated  3p^  'S)  [10].  However  this  state  has  odd  parity,  opposite 
to  the  ground  state’s  even  parity.  Therefore,  if  we  consider  only  dipole  transitions,  we  must 
consider  a  three  photon  process.  However  a  3.9  eV  photon  does  not  propagate  long  distance  in 
the  atmosphere,  so  we  do  not  consider  this  process. 

The  next  even  state  is  designated  4p[l/2]  and  has  an  energy  of  13  eV,  so  it  can  be  excited  by  a 
four  photon  process  with  a  3.25  eV  photon.  Then  the  excited  state  can  be  ionized  in  a  single 
photon  process,  which,  as  we  will  see,  goes  very  quickly.  Hence  the  4  photon  excitation  of  argon 
is  the  rate  limiting  step.  Since  it  is  a  resonant  process,  and  the  competing  process,  four  photon 
direct  ionization  of  O2  is  not,  we  would  expect  it  to  go  faster.  Thus  one  could  envision  a 
cascade  process.  The  most  difficult  first  task  is  to  excite  the  argon  atom  with  4  photons.  Then 
ionizing  it  with  a  single  photon  is  much  easier. 

Although  we  focus  on  the  four  photon  excitation  of  the  4p[l/2]  state  with  a  3.25  eV  photon,  there 
are  other  possibilities  as  well.  The  4s [1 1/2]°  state  could  be  excited  by  5  photon  process  with  a 
2.33  eV  photon,  followed  by  two  photon  photo  ionization  of  the  excited  state.  Also  this  state 
could  be  excited  by  a  four  photon  process  if  one  of  the  intermediate  state  matrix  elements 
happens  to  be  unexpectedly  large  for  a  quadrupole  transition.  If  this  occurs,  it  turns  out  that  there 
are  several  other  5p  states  just  2.9  eV  above,  so  it  could  be  excited  by  a  single  phonon,  followed 
by  single  photon  excitation. 

To  quantify  this  theoretically  however,  one  would  have  to  do  multi  photon  quantum  mechanical 
perturbation  calculations  for  many  electron  atoms  and  molecules.  We  do  not  attempt  that  here. 
Alternatively,  one  could  attempt  to  settle  it  experimentally.  Using  a  tunable  high  power  short 
pulse  laser  (EEL?  dye  laser?,  frequency  doubled  Ti:saphire  laser?),  one  could  attempt  to  ionize 
air  with  photon  energies  around  3.25  eV  and  investigate  whether  there  is  a  resonant  process  to  be 
exploited.  Furthermore  one  could  also  perform  analogous  experimental  investigations  for  photn 
energies  around  2.33  and  2.9  eV. 
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While  the  necessary  theoretical  calculations  are  difficult,  one  can  obtain  considerable  insight  by 
doing  a  variety  of  simpler  quantum  mechanical  calculations.  In  fact,  these  simple  calculations  do 
indeed  indicate  that  ionization  via  the  argon  channel  is  several  orders  of  magnitude  faster  than 
the  direct  4  photon  ionization  of  O2.  In  Section  II  we  look  at  two  state  systems,  in  Section  m, 
multistate  systems  (Fermi’s  golden  rule).  In  Section  FV  we  look  into  single  photon  excitations  of 
atoms  and  molecules  and  give  rates  for  simple  systems  or  rough  approximations  for  more 
complicated  ones.  In  Section  V  we  look  into  single  photon  ionization  ionization,  while  in 
Section  VI  we  look  into  multiphoton  process  for  atoms  and  molecules.  Section  Vn  discusses 
solutions  to  the  rate  equations  for  4  photon  ionization  of  argon,  and  Section  VUI  summarizes. 


II.  The  Two  Level  System: 


This  section  and  the  next  briefly  reviews  quantum  mechanical  perturbation  theory.  These  are 
rather  brief,  as  this  information  can  be  found  in  the  standard  texts.  The  system  is  governed  by 
Schroedinger’s  equation 


at 


(1) 


where  H  is  the  Hamiltonian  operator.  If  H  were  time  independent,  we  would  have 


\j/(r,t)  =  i|/(r)exp- 


(2) 


Let  us  say  that  in  the  system  we  consider,  there  are  only  two  states,  denoted  by  subscripts  zero 
and  one.  They  are  orthogonal  to  each  other  and  each  is  normalized  to  unity  in  the  usual  way. 
since  the  system  is  a  two  state  system,  these  eigenfunctions  are  also  complete.  Then  the  wave 
function  will  always  be  the  linear  combination 


\|/(r,t)  =  Co\|/o(r)exp- 


ri 


iE„ 


/: 


t  +  C,\|/,(r)exp- 


iE, 


(3) 


where  because  of  the  completemess,  |Co^|  +  |Ci^|  =  1.  Let  us  say  now  that  there  is  a  time 
dependent  perturbed  Hamiltonian,  but  this  perturbation  does  not  change  the  two  state  nature  of 
the  system.  Then  the  solution  to  Schroedinger’s  is  as  written  in  Eq.(3),  except  that  now,  each  C 
must  be  a  function  of  time.  Assuming  this,  and  further  assuming  that  the  time  dependence  of  the 
perturbation  goes  as  EoCos  tot,  we  substitute  into  Schroedinger’s  equation.  On  the  left  hand  side 
we  find  terms  proportional  to  the  >j/’s  and  their  time  derivatives.  Also  many  terms  cancel 
because  the  tj/’s  are  solutions  to  the  time  independent  Schroedinger’s  equation.  Multiplying  the 
resulting  equation  by  first  \j/o  and  integrating  over  space,  and  then  \|/],  we  find  two  equations  for 
the  time  derivatives  of  the  C’s.  These  are 
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(4a) 


dQ 

— -  =  — C,  exp-  lAcot 
dt  \h 

dC,  .... 

— ^  = - Cfl  exp  1  Aa)t  (4b) 

dt  ih 

where  Mco  =  E,  -E^  -to  and  M=<:\(/o*  |6H\|/i>=/d\\|/o*  5H\|ri.  In  writing  Eqs.(4a  and  b)  we 

have  used  the  Hermitian  properties  of  5H.  Also  we  assume  that  Ao)  is  small  (leading  to  secular, 
long  time  scale  changes  in  the  C’s)  and  have  thereby  also  dropped  the  higher  frequency  term  on 
the  right  and  side.  Defining  =  |M^|/  ^  Eqs.(4)  are  harmonic  oscillator  equations  with  a 
frequency  Q  given  by 

(5) 


If  Ao)  =  0,  the  perturbation  causes  the  system  to  oscillate  between  states  0  and  1  with  frequency 

Q. 


III.  The  Multistate  System,  Fermi’s  Golden  Rule: 

Let  us  imagine  the  lower  state,  zero  still  being  a  single  state,  but  the  upper  state  being  one  of 
many  possible  states.  That  is  we  replace  the  subscript  1  with  j,  which  is  now  an  index  of 
summation.  Equation  (4b)  can  be  integrated  once  to  give 


iMi*C„  .  .sm(Aa).t/2) 

C:  = — — ^expi(Aa)it/2) — 

'  h  y  \  i  )  Aa).t/2 


(6) 


Multiplying  the  equation  for  Cj  by  the  equation  for  dCj*/dt  and  adding  the  complex  conjugate,  we 
find 
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^  =  4nQ.;6(Aa.,)|C.r  (7) 

where  we  have  used  the  convention  that  for  long  time,  a)‘'*sina)t-^  7t6(a)).  Of  course  because  of 
the  completeness,  we  still  have  |Cof  +  2j  |Cjf  =  1.  Hence  to  find  the  loss  from  the  original  state, 
we  sum  over  all  j  and  find 

^l^  =  -24n|a.iC,f8(Ato,)  (8) 

which  is  valid  as  long  as  |Cof  »  |Cip.  The  delta  function  only  makes  sense  for  an  integral,  not 
a  summation.  To  convert  the  summation  to  an  integral,  we  need  a  density  of  states  p(o)). 
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(9a) 


^  =|da)p({D) 


and  this  depends  on  the  details  of  the  physical  situation.  Integrating  we  find 


dl 


-SHn^c,fp(A(Oj=o)- 


-471 


A(On 


(9b) 


where  A(Ob  is  the  frequency  width  of  the  band  of  states.  We  will  denote  this  a  molecular,  as 
opposed  to  an  atomic  process,  for  reasons  which  will  become  clear  shortly. 

Alternatively,  if  the  state  is  isolated  but  the  laser  has  a  spread  in  frequency  5(0,  so  that  the 
irradiance  I  =  JdtoI((D),  then  the  sum  over  j  in  Eq.  (8)  becomes  an  integral  over  O).  If  the 
irradiance  has  frequency  spread  8co,  which  is  many  instances  is  the  reciprocal  of  the  pulse  time 
then  the  rate  goes  roughly  as 


M 

dt 


-471Q„  Tp 


(10) 


We  call  this  the  rate  for  an  atomic  process  since  these  states  are  isolated. 

If  the  Q,^  for  the  atomic  and  molecular  processes  are  about  the  same,  to  optimize  the 
irradiance  for  maximum  excitation  during  the  pulse,  we  set  47ti2jp  ~  1  (larger  irradiance  will 
cause  the  excitations  to  oscillate  during  the  pulse).  At  this  irraadiance,  the  molecular  excitation 
rate  is  smaller  by  a  rate  of  47tA(0„x  ,  or  at  least  three  orders  of  magnitude  for  a  1  eV  band  and  a 
100  fs  pulse. 


IV.  Single  Photon  Excitation  for  Atoms  and  Molecules: 

Here  we  calculate  for  simple  systems,  and  estimate  for  more  complicated  systems,  the  single 
photon  excitation  rate.  The  perturbed  Hamiltonian  is 

6H(r,  t)  =  -er»e„  cos  tot  (11) 

that  is  the  energy  of  the  displaced  charge  in  the  oscillating  field.  The  radiation  field  has  a  spatial 
dependence  as  well,  but  the  wavelength  is  normally  large  compared  to  atomic  or  molecular 
distances  so  we  neglect  it  in  the  dipole  approximation.  Since  in  the  derivation  in  Sections  2  and 
3,  the  perturbed  time  dependence  was  expressed  in  term  exponentials  rather  than  cosines,  the 
matrix  element  is  <ai/i*|er»e\|/o>  will  have  an  additional  factor  of  */2. 

IVa.  Atoms:  For  the  two  states,  we  use  the  ground  state  and  1P°  excited  state  for  the  hydrogen 
atom.  These  wave  functions  are 
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(12a) 


¥o  = 


J 


V,  = 


lcos0 


J 


(12b) 


where  0  is  the  polar  angle  to  the  z  axis,  defined  by  the  direction  of  the  electric  field.  The  quantity 
ao  is  the  Bohr  radius  ao  =  ^  ^/me^.  The  integrals  are  all  straightforward,  and  we  find  the  matrix 
element  is  given  by  M  =  l.beEoao,  and  the  oscillation  frequency  between  states  by 


n 


(13) 


In  cgs  units,  take  the  approximate  values  e  =  5xl0  ’®,  ^=10’^^,  ao  =  5x10'^.  For  the  electric  field 

10  0  4  10  1  /o 

we  assume  a  canonical  irradiance  of  10  W/cm  .  In  this  case,  e  =  6x10  (I/IO  )  ,  so  O  » 
2.4xl0''‘(F10‘^)''^. 


rVb.  Molecules:  For  molecules,  the  excited  states  are  not  isolated  as  they  are  in  an  atom.  Rather 
they  occur  in  bands  because  for  each  electronic  excitation,  there  are  a  very  large  number  of 
vibrational  and  rotational  state  states  as  well.  Schematically,  the  electronic  excitations  are 
sketched  below  as: 
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An  actual  measured  spectrum  for  the  nitrogen  molecule  between  3000  and  4800  angstroms  is 
shown  below  as  well. 


Clearly  there  is  no  single  excited  state,  but  rather  nearly  a  continuum  of  states  occurring  in  bands. 
Furthermore,  for  the  low  level  states,  the  bands  do  not  overlap,  so  it  is  not  possible  to  excite  and 
ionize  with  a  cascade  of  low  energy  photons.  If  the  band  has  frequency  width  Acob,  and  if  the 
matrix  elements  for  all  transitions  within  the  band  are  about  the  same,  we  can  consider  the  band 
to  be  a  macro  state  with  frequency  width  Ao)b,  meaning  a  density  of  states  AcOb'*  if  the 
summation  over  states  is  converted  to  an  integral  over  frequency.  Hence  the  transition  rate  as 
given  in  Eq.  (9b) 

This  idea,  that  the  molecular  excitation  rates  are  expected  as  much  as  three  orders  of  magnitude 
smaller  than  resonant  atomic  rates  is  supported  by  the  simple  calculations  done  here.  It  is  the 
basis  for  our  expectation  that  photo  ionization  of  the  atmosphere  may  proceed  more  easily 
through  resonantly  enhanced  excitation  in  the  argon  atom,  than  from  analogous  processes  in  the 
oxygen  molecule. 
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V.  Single  Photon  Ionization: 

Here  we  do  not  distinguish  between  ionization  of  atoms  and  molecules,  because  in  each  case,  the 
final  state  is  similar,  a  free  electron.  If  the  free  electron  has  high  energy  compared  to  the 
ionization  energy,  it  is  not  strongly  affected  by  the  nucleus.  We  assume  that  in  this  case,  the 
electron  is  a  freely  propagating  wave.  If  it  is  in  a  box  of  size  L  on  a  side,  the  wave  function  for 
the  free  state  is  expik»r  in  the  high  energy  limit.  Then  the  matrix  element  is  given  by 

M(k)  =  — ■j=r^ - sTtIO*  Texp  [ik»r  -  r  /  a„  ] + c.c.  (14) 

2V7t(a„L) 

where  again,  Eq  is  the  amplitude  of  the  electric  field  expressed  as  a  cosine.  Now  say  0  is  the 
angle  between  e  and  k,  0”  is  the  angle  between  r  and  e,  and  0’  is  the  angle  between  r  and  k.  If  we 
let  k  be  the  polar  axis,  then  we  can  relate  0”  to  the  other  angles  cos0”  =  cos0’cos0  + 
sin0’sin0cos((l)’-x)  where  (j)’  and  %  are  the  azimuthal  angles.  The  volume  element  is  then 
r^drsin0’d0’d(t).  The  integrals  over  (|)’  vanish,  so  we  are  left  with 

M(k)  =  p— - — jp  f  f  r^dr  sin  0 '  d0 '  ee^r  cos  0  cos  0  'exp  [ikr  cos  0'-r/ap]+ c.c. 

V7i(a„L) 


_  167te„ka/cos0 

(l  +  kV)' 

To  get  the  ionization  rate,  we  use  the  golden  rule,  Eq  (8).  To  so  so  we  need  the  density  of  states. 
The  energy  of  an  electron  with  wave  number  k  is  Ek  k^/2m,  so  the  density  of  states  with 
energy  Ek  is 


P(E0  = 


lV^^VeT 


dE, 


Using  the  relation  between  energy  and  frequency,  we  can  do  the  integral  in  Eq.  (8)  to  get  the 
ionization  rate  w. 


256mk^e^e/a  ’ 
w  = - V 

3^^(l  +  kX^) 


(15) 


where  we  have  integrated  over  angle.  Of  course  k  is  related  to  the  photon  energy  and  ground 
state  energy  Ep  (a  negative  number  by)  /  2m  =  ftto+ E„  .  Asa  function  of  k,  w  in  Eq.  (15) 

has  a  maximum  at  kap  =  1,  at  which  point,  we  find  w  ~  10*^  s'^  for  our  canonical  irradiance  of 
lO’^  W/cm^.  For  larger  k,  the  rate  falls  off  rapidly  with  photon  energy. 
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Equation  15  also  shows  that  at  threshold,  k  — >  0,  the  rate  rapidly  goes  to  zero.  This  is  not  the 
correct  behavior,  and  it  arises  from  the  assumption  of  a  freely  propagating  electron  wave.  Near 
zero  energy,  the  electron  wave  function  is  strongly  affected  by  the  nucleus.  The  calculations  of 
ionization  rate  near  the  absorption  edge  are  difficult  [12].  To  get  an  approximate  value,  we 
assume  that  the  matrix  element  and  state  density  do  not  change  very  much  on  going  just  below  to 
just  above  the  ionization  energy.  This  is  rather  like  the  orbits  of  particles  near  the  separatrix  in  a 
sinusoidal  potential.  Other  than  the  fact  that  the  just  barely  trapped  orbits  are  confined  and  the 
just  barely  untrapped  orbits  are  not,  the  nature  of  the  orbits  are  very  similar  along  any  nearly 
portion  of  the  orbit. 

Hence  we  approximate  the  ionization  rate  with  the  transition  rate  to  very  high  Rydberg  states  in 
the  atom.  Since  the  ground  state  has  1=0,  the  only  allowed  dipole  transitions  are  to  1=1  states.  In 
terms  of  the  ionization  energy  Ei  =  hdi- ,  the  energy  of  a  Rydberg  state  of  quantum  n  is  -EJn^.  If 
the  radial  part  of  the  wave  function  is  R(r),  Schroedinger’s  equation  for  R  becomes 


1  d  2dR 
r^dr^  dr 


(16) 


where  we  assume  n  is  large.  For  large  r,  the  angular  momentum  barrier  (last  term  on  rhs)  plays 
almost  no  role,  and  for  small  r,  the  energy  eigenvalue  (the  first  term  on  rhs)  plays  almost  no  role. 
If  one  uses  a  WKB  approximation  for  R,  the  function  is  oscilliatory  in  r  as  long  as  ao  <  r  <  n^ao, 
and  it  exponentially  damps  away  for  larger  r  (the  eigenvalue  is  selected  to  eliminate  the 
exponentially  growing  solution).  Thus  for  r<n^ao,  the  magnitude  of  R  is  roughly  uniform,  so  that 
the  normalization  for  the  wave  function  is  roughly  R  ~  l/(ao^^^n^).  The  value  of  the  radial  part  of 
the  ground  state  is  given  in  Eq.  (10a).  so  the  matrix  element  is  given  roughly  by  M  ~ 

The  ionization  rate  is  given  by  Eq.  (8),  which  is  roughly 

w^47i(n.fe-^-  (17) 

Ei 


where  we  have  used  for  the  density  of  states  for  high  Rydberg  states 

dn((o)  _  1  _  n^ 

do)  2  V  (0^  2(0; 

If  the  initial  state  were  not  the  ground  state,  but  a  higher  Rydberg  state  with  principle  quantum 
number  no,  an  analogous  argument  shows  that  the  matrix  element  is  multiplied  by  and  the  rate 
goes  roughly  as 


w  «  4n(i),nJ‘ 


(17a) 


Notice  that  at  threshold,  the  rate  no  longer  vanishes,  as  the  high  energy  approximation,  Eq.  (15) 
would  have  predicted.  Instead,  at  threshold,  one  has  an  ionization  rate  given  roughly  by  Eq.  (17). 
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As  the  photon  energy  increases,  the  rate  decreases,  ultimately  being  given  by  the  high  energy 
approximation,  Eq.  (15).  Expressing  the  rate  in  terms  of  an  ionization  cross  section,  w  = 
ai[I/ So);  ],  and  expressing  the  irradiance  as  I  =  cEo^/Sti,  we  find  near  the  ionization  threshold. 


a,  =32Tt' 

V' 

2 

or 

G:  ~32n^nJ‘ 

0 

the  latter  being  for  an  excited  state.  For  the  ground  state,  we  find  Oi  ~  5x10  '^  cm^,  about  the 
right  value.  Thus  approximating  the  wave  function  for  energy  just  above  zero,  with  that  for 
energy  just  below  zero  gives  a  simple  and  reasonably  accurate  estimate  of  the  ionization  cross 
section,  or  ionization  rate,  just  above  threshold. 


VI.  Multi  Photon  Processes: 

Multi  photon  processes  are  extremely  complicated  to  calculate,  even  for  a  hydrogen  atom,  so  we 
will  confine  ourselves  to  their  rather  general  properties.  We  find  that  if  the  levels  are  isolated,  as 
in  an  atom,  there  are  still  oscillations  between  the  states.  While  the  process  is  nonlinear  in  the 
matrix  element  (i.e.  in  e),  it  is  still  linear  in  the  two  states,  zero  and  one,  which  are  connected  by 
the  multi  photon  process.  On  the  other  hand,  if  the  final  state  is  a  band  of  states  (or  ionization), 
the  basic  behavior  is  again  as  in  the  case  of  single  photon  excitation.  That  is  there  is  an  analog  to 
the  golden  rule,  but  the  coefficient  how  goes  as  some  higher,  but  even  power  of  8.  For  multi 
photon  processes,  it  now  becomes  essential  to  consider  other  excited  states.  We  will  define  by 
zero  and  one  the  two  isolated  states  we  expect  to  connect  with  a  multi  photon  process  and 
consider  other  states  as  well  with  an  index  j  which  goes  from  2  to  some  upper  limit  P.  The  P+2 
eigenfunctions  once  more  are  assumed  to  form  a  complete  orthonormal  set.  Then  Eqs  4  become 

^  =  ^C,  exp-iAa)<„t  +  exp-iAa)„jt  (18a) 

dt  "  i«  ' 

^  ^  Co  exp-  i  A(0,„t  +  ^  C j  exp-  i Aco, ^t  (1 8b) 

dt  ^  ih  ‘ 

dC:  M.  Mi,  P  M,,^ 

— -  =  — exp-iAo),  t  + — Ac,  exp-iAa)j,t  +  V — exp-iAtOji^t  (18b) 
dt  ift  \h  ^  Hj  i^ 

where  now  the  Ato’s  are  defined  analogously,  except  that  now  A(0jk  =  -AtOicj  and  the  matrices  have 
the  same  symmetry  properties  as  before.  To  illustrate  the  nature  of  the  coupling,  let  us  assume 
that  states  0  and  1  have  an  energy  separation  close  to  2  ^O) ,  and  assume  further  that  only  a  single 
other  state,  say  Cj  is  important  for  the  coupling.  Then  integrating  Eq.  (18b)  we  find 

M  M 

Cj  =  :  - -  Co  exp-iAco^^t  exp-iAm,,!  (19) 

^AtOjo  SAo)2, 
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where  in  Eqs  (18  and  19)  the  Ato’s  are  not  small.  Then  we  substitute  for  C2  on  the  right  hand 
sides  of  Eqs  (18a  and  b)  to  determine  the  next  order  equation  for  Co  and  Cj .  Then  we  get  for  Ci 


“IT  "  -J  a  '  ^0  exp-  i(Aco,2  +  Ato^o  )t 

dt  AtOjo 


(20) 


and  analogously  for  Cq.  In  Eq.  (20)  we  have  kept  only  those  exponentials  which  are  at  low 
frequency,  that  is  at  a  frequency  near  [(Ei  -  E2)/  ti  -  2(i)]  =  0.  Equation  20  and  the  analogous 
equation  for  Co  are  still  linear  in  the  C’s  even  though  the  coupling  coefficient  is  now  proportional 
to  rather  than  simply  £.  They  are  still  harmonic  oscillator  equation  equations,  so  the  process 
still  has  the  populations  oscillate  sinusoidally,  but  now  at  frequency  Q2  =  |Mi2M2o/^  ^Aa)|, 
between  the  two  states  which  are  connected  by  the  two  photon  process.  The  rate  for  the  two 
photon  process  as  compared  to  the  one  photon  process  is  given  roughly  by 


onephotonrate  0-2  M  eea„ 

twophotonrate  Q,  /?A(0  fiAo) 

This  ratio  is  proportional  to  the  electric  field  e,  and  the  rate  itself  is  proportional  to  8  . 

Two  photon  processes  (and  even  more  emphatically,  greater  than  two  photon  processes)  are  not 
as  constrained  by  selection  rules  as  are  single  photon  processes.  For  instance,  states  zero  and  one 
may  be  forbidden  dipole  transitions.  However  each  may  have  a  non  zero  matrix  element  between 
itself  and  state  two,  and  there  are  many  possible  intermediate  states.  Thus  usually  we  would 
expect  a  multiphoton  path  between  the  two  states.  For  more  than  two  photon  processes,  there  are 
even  more  channels.  Nevertheless,  if  only  dipole  transitions  are  considered,  odd  to  even  states 
must  be  connected  by  odd  numbers  of  photons,  and  visa  versa. 

If  the  intermediate  state  is  a  band  of  states  rather  than  an  isolated  state,  the  transition  rate  will  not 
be  significantly  different,  the  distinction  only  arises  if  the  band  is  resonant  or  nearly  resonant. 
However  since  Ao)o2  is  not  small,  the  intermediate  state  can  be  an  isolated  state  or  a  band  of 
states.  But  if  the  final  state  is  a  band  of  states,  the  result  is  still  significantly  different.  Say  that 
states  zero  and  two  are  isolated  (but  two  may  be  part  of  a  band  without  significant  effect)  and  one 
is  a  band  of  nearby  states.  In  that  case  the  procedure  is  the  same  as  in  Section  4b.  We  solve  for 
Cl  by  integrating  Eq(20)  once  in  time.  Multiply  this  by  the  complex  conjugate  of  Eq.  (20),  and 
add  the  complex  conjugate  of  the  product  to  get  an  equation  for  the  rate  of  change  of  |Cip.  It  is 


=  471 


/i^AtO^o' 


(22) 


Thus  if  the  final  state  is  a  band  of  states,  of  bandwidth  AcDb,  and  all  these  states  have  about  the 
same  matrix  element,  we  may,  as  for  the  single  photon  case,  consider  it  to  be  a  single  macrostate 
in  a  frequency  range  AcOb,  giving  a  density  of  states  A(0b*\  so  that 
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dt  /^''AcjOjo^AcOb 


Q. 

4n - ^ 


=  471 


AcD^AWg  A(Og 


(23) 


Instead,  if  the  states  are  isolated  but  the  laser  has  bandwidth  6(Ol,  possibly  equal  to  Xp'’,then 


d|C, 

dt 


■~4nQ.2  \ 


(24) 


Once  again,  if  Q,2'^p  ^  1,  the  incoherent  rate  is  not  that  much  less  than  the  coherent  rate.  If  we 

define  Qn  =  (J0L[eeao/E]",  where  E  is  a  typical  energy  gap  between  the  relevant  states,  then  the  n 
photon  atomic  excitation  rate  goes  roughly  as  47ii2„^Xp.  Similarly,  the  n  photon  ionization  rate 
from  the  ground  state,  by  analogy  to  Eq.  (17),  goes  roughly  as 

w„  ~  47tO)L[eEao/E]^".  (25) 

Thus  the  n  photon  excitation  rate  exceeds  the  n  photon  ionization  rate  by  a  factor  of  47ia)LTp,  a 
factor  of  several  orders  of  magnitude.  This  then  is  the  principle  motivation  for  our  expectation 
that  four  photon  photo  excitation  of  the  argon  might  be  a  much  faster  process  than  direct  four 
photon  ionization  of  O2.  The  former  is  a  resonant  process,  the  latter  is  not. 


VIL  The  Rate  Equations: 

We  would  like  to  finish  by  writing  out  rate  equations  for  the  ground  state,  excited  state  and 
electrons.  Using  subscripts  g,  x  and  e,  we  find 

^  =  ^x(n*-ng)  (26a) 

^  =  -^x(nx-ng)-Vin,  (26b) 


dt 


(26c) 


where  is  the  four  photon  excitation  rate  for  argon,  and  Vj  is  the  ionization  rate  of  the  excited 
argon  atom.  Also  we  have  included  the  fact  that  the  electric  field  excites  the  atom  and  de  excites 
it  at  the  same  rate.  Equations  give  the  trivial  conservation  law  that  total  number  of  particles  is 
conserved.  We  expect  that  Vi »  so  that  the  excited  state  is  quasi  steady  in  time.  Setting  the 
left  hand  side  of  Eq.  (26)  to  zero,  we  can  solve  nx  =  [i2x/(  Vi+Qx)]ng  «  ng.  Hence  the  ionization 
rate  is  simply  the  excitation  rate  of  the  argon.  Because  it  is  a  resonant  process,  we  expect  this 
rate  to  be  much  larger  than  the  analogous  ionization  or  excitation  rates  for  the  molecular 
components  of  the  atmosphere. 
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VIQ.  Summary. 

In  this  memo,  we  have  laid  out  the  case  that  the  optimum  way  to  ionize  the  atmosphere,  at  a 
controlled  range,  with  a  short  pulse,  high  power  laser  may  well  be  to  ionize  the  argon  rather  than 
the  oxygen  or  nitrogen.  The  reason  is  that  the  argon  is  an  atom,  not  a  molecule.  As  such  its 
electronic  excitations  are  isolated  rather  than  part  of  a  band  structure.  If  the  laser  is  properly 
tuned,  they  should  be  much  easier  to  excite.  We  have  also  laid  out  a  possible  route,  a  4  photon 
excitation,  followed  by  photo  ionization  of  the  excited  atom.  Other  routes  are  also  possible.  Our 
rough  estimate  is  that  this  rate  could  be  greater  than  the  direct  four  photon  ionization  of  oxygen 
by  as  much  as  several  orders  of  magnitude.  To  evaluate  this  more  accurately,  we  would  have  to 
do  the  multi  photon  multi  electron  calculations  not  only  for  argon,  but  also  for  the  oxygen 
molecule.  That  way,  we  could  compare  the  two  and  see  which  is  better. 

An  important  issue  is  that  even  if  we  do  these  calculations,  would  we  believe  them?  To  this 
author  they  look  so  complicated  that  I  would  not  be  so  sure.  Of  course  if  we  do  adopt  such  an 
ionization  scheme,  we  would  ultimately  have  to  do  them.  However  for  the  short  or  intermediate 
term,  an  alternative  could  be  to  simply  do  the  experiment. 

One  possibility  would  be  to  use  our  own  Tiisaphire  laser  and  frequency  double  it.  Alternatively, 
at  the  Jefferson  Lab  FEL,  with  frequency  multiplication,  they  could  be  in  the  optical  regime, 
albeit  at  rather  low  power  and  efficiency  [13].  However,  within  the  next  six  months  to  year,  they 
do  expect  to  have  their  system  operating  much  more  efficiently  in  the  uv  and  optical  range. 

Either  way,  the  idea  would  be  to  the  send  the  laser  light  into  a  target  chamber  containing  only  air, 
but  possibly  at  reduced  pressure,  so  that  if  there  is  ionization,  the  electron  lifetime  would  be 
longer  and  easier  to  measure.  Then  we  would  look  at  the  ionization  as  a  function  of  frequency 
and  irradiance.  The  irradiance  one  could  probably  most  easily  vary  by  changing  the  focal 
conditions  of  the  beam.  We  would  like  to  vary  the  photon  energy  between  about  3  and  3.5  eV, 
around  the  upper  limit  of  photon  energy  that  propagates  in  the  atmosphere,  and  see  of  there  is  a 
resonant  process  at  a  photon  energy  of  about  3.25  eV.  If  so,  it  would  indicate  that  four  photon 
excitation  of  argon  is  playing  the  key  role  in  the  ionization  process.  Other  photon  energies  which 
might  be  of  interest  for  ionization  of  the  argon  could  be  2.33  and  2.9  eV. 
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